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Abstract 

It is shown that the exchange-correlation part of the action functional 
Axcipi^^t)] in time-dependent density functional theory , where p{r,t) is the 
time-dependent density, is invariant under the transformation to an accel- 
erated frame of reference p{r,t) — > p'{f,t) = p{f+x{t),t), where x{t) is 
an arbitrary function of time. This invariance implies that the exchange- 
correlation potential in the Kohn-Sham equation transforms in the following 
manner: Vxc[p'',r,t] = Vxc[p',r + x{t),t]. Some of the approximate formulas 
that have been proposed for Vxc satisfy this exact transformation property, 
others do not. Those which transform in the correct manner automatically 
satisfy the "harmonic potential theorem" , i.e. the separation of the center of 
mass motion for a system of interacting particles in the presence of a harmonic 
external potential. A general method to generate functionals which possess 
the correct symmetry is proposed. 
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Time dependent density functional theory (TDFT) |T|-§| is a very valuable tool for the 
study of the dynamic behavior of correlated electronic systems under the influence of external 
probes which can be represented as time-dependent potentials. The range of problems to 
which this formalism is successfully applied is expanding 0]: it includes the calculation 
of atomic collision cross sections, the behavior of atoms under very strong electromagnetic 
fields, the calculation of excitation energies in atoms, molecules, and solids . Crucial to 
the success of the theory is the availability of simple approximations for the exchange- 
correlation potential Vxc[p', r, t], which is an ordinary function of r and t, but a complicated, 
non-local functional of the density distribution p(r^, t') for times t' earlier than t. Examples 
of approximations for V^c are the "adiabatic time dependent local density approximation" 
(ATLDA) of Zangwill and Soven and the local linear response approximation, with a 
frequency-dependent exchange-correlation kernel, of Gross and Kohn p[. In constructing 
such approximations, one should try to make sure that the approximate form satisfies as 
many of the known exact properties of the functional as possible. One example of rigorous 
constraints on admissible approximate functionals is provided by the scaling relations for 
the static density functional theory Another type of constraint, which is peculiar to 
the time-dependent theory, has been recently pointed out by Dobson [0. He has observed 
that, according to a general "harmonic potential theorem" (HPT), a system of interacting 
electrons confined by an external harmonic potential Vo(r) = (l/2)r ■ K ■ r, must possess 
solutions in which the dynamics of the center of mass is completely decoupled from that of 
the internal degrees of freedom. Assuming that the system is initially in the ground state, 
the center of mass wave function is a minimum uncertainty gaussian centered about the 
origin. When a time-dependent uniform electric field E(t) is applied to such a system, the 
time dependence of the center of mass wave function is such that the gaussian is rigidly 
transported (to within a time-dependent phase factor), the position of its center being 
determined by the solution of the classical equation of motion 
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where B is a. uniform and constant magnetic field. Tlie fact that the "relative" dynamics 
is not affected by the external force implies that the electronic density must also be rigidly 
transported, following the motion of the center of mass, i.e., 

p{r,t)=po{f-RcMit)) (2) 

where poir) is the static electronic density in the absence of the driving field. Of course, 
one would hope that approximate versions of TDFT were able to reproduce this elegant 
exact result. Dobson's discovery ^ is that this is not always the case. For example, the 
Zangwill-Soven adiabatic ATLDA does satisfy this requirement, but the Gross-Kohn 
linearized response theory, which includes a frequency dependence (i.e., a memory) in 
the exchange-correlation kernel, does not. Other approximations, based on inhomogeneous 
hydrodynamics, also fail to satisfy the HPT . The root of the difficulty was traced by 
Dobson to the inclusion of non-locality in time but not in space. He developed a simple 
modification of the TDLDA, whereby the density modulation induced by the driving field is 
decomposed in two parts, one corresponding to a rigid shift of the whole distribution (i.e, a 
motion of the center of mass), the other including the remaining internal deformations. The 
Gross-Kohn formalism is applied only to the latter part, while the center of mass component 
is still treated in the adiabatic LDA. This formalism satisfies the HPT. 

In this paper, we demonstrate that the occurrence of the HPT in TDFT is in fact the 
consequence of a simple symmetry of the exchange-correlation part of the action. The 
approximations that Dobson finds to satisfy HPT are the ones that satisfy this symmetry, 
and those that violate the HPT are the ones that violate this symmetry. Furthermore, the 
cure he proposes to make the TDLDA consistent with the HPT can be viewed as an instance 
of a general procedure whereby the correct symmetry of the exchange-correlation potential 
can be built at the outset in any approximation. 

We now present our results in detail. The basic theorem of TDFT is the Runge-Gross 
theorem, according to which, for any time-dependent density p(r, t) defined within an ap- 
propriate class of functions in the interval to < t < ti, there exists an essentially unique time 
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dependent potential V{f, t), such that the solution \ip{t)) of the time-dependent Schroedinger 
equation 

{t^^-Ho-V{t)mt))=0, (3) 

with the condition that ip(to) is the Hohenberg-Kohn ground-state wave function correspond- 
ing to the density po{r) = p{f,to), yields p{f,t) as the expectation value {ip{t)\p{'r)\ip{t)) of 
the density operator p{f) — J2iLi S{f— fj). Here 

1=1 * J ' ' 

is the familiar Hamiltonian for N interacting electrons in the presence of a uniform magnetic 
field B described by the vector potential A{r) — B x f/2, and 

v{t)=j:vin,t) (5) 

i=l 

is the time-dependent potential. The expression "essentially unique" refers to the possibility 
of adding to V{f, t) an arbitrary function f(t) of time, which causes the wave function to be 
multiplied by a phase factor exp[i J f{t)dt] without affecting the density. The "appropriate 
class" of densities includes all the densities for which the potential V{f, t) exists, and is Taylor 
expandable, with finite convergence radius, in a neighborhood of to. Under the assumption 
that at time to the system is in the ground state - an assumption that we make in the rest 
of this paper - the Runge-Gross theorem establishes the existence of an essentially unique 
mapping from time-dependent densities to time-dependent potentials. 

Let us now look at the density p{f, t) from the point of view of an accelerated observer 
whose position, relative to the original reference frame, is given by the function x{t). It 
is assumed that the accelerated observer uses cartesian axes that remain parallel to the 
axes of the original reference frame, i.e., there is no rotation. The density seen by the 
observer, in terms of his own r coordinate is p'{r,t) — p{f + x{t),t). It is assumed that 
x{to) = dx{t)/dt{t = to) = so that both the density and the wave function seen by 
the accelerated observer coincide with those seen by the inertial observer at time to- Our 
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first important result is that there exists a time-dependent potential V^'(r, t) that generates 
p'{r,t) when the Schroedinger equation is solved. The explicit form for V'{f,t), aside from 
an inessential additive function of time, is 

N N 

V\r,t) = Vif+ f(t),t) + md{t) ■ 5]r, + {e/c){v{t) x B) ■ ^rl, (6) 

1=1 i=l 

where a{t) and v{t) are the second and first derivatives of x{t) with respect to time. The 
wavefunction ip' corresponding to this transformed potential is related to the wavefunction 
ip corresponding to potential V by the transformation 

\ij'{t)) = uiAit)\m) (7) 

where the unitary transformation Ui{t) is defined as 

c -> 

Ui{t) = exp[—ifi ■ {mv{t) — —B x x{t))]exp[ipi ■ x{t)]. (8) 

In writing the above formulas we have assumed that there is a uniform and constant 
magnetic field B acting on the electrons. This is because we want to discuss later the 
generalized harmonic potential theorem in the presence of a magnetic field. The familiar 
situation is recovered by putting B = in the above formulas. The detailed mathematical 
proof of the key equations (|^) and (0) is provided in the appendix. However, the correctness 
of these equations is physically evident. Equation @ says that the potential seen by the 
accelerated observer is the original potential expressed in terms of the new coordinate, plus 
an inertial force —ma{t) that couples like a uniform gravitational field to the center of mass 
of the system, plus a uniform electric field v x B/c that arises from the transformation 
of the uniform magnetic field. Equation says that the new wavefunction is obtained 
from the former by translating the coordinate of each electron by x{t) and the momentum 
of each electron by mv{t). When the magnetic field is present an additional translation of 
the momentum by —eB x x{t)/2c is required in order to restore the original gauge for the 
vector potential. Equations and (J^) enable us to determine the transformation of the 
"internal" action functional, defined as follows: 

5 



Mp]= j\Hm^^-Ho\m)dt (9) 
(Note the absence of the external potential in this definition). We obtain 
A[p'] = A[p] + N f\ma{t) ■ Rcuit) + -{m x B) ■ Rcuit) + + ^{B x x{t)) ■ v{t)]dt, 

Jto C Z Zc 

(10) 

where 

NRcM{t) = I fp{f,t)df (11) 

is the coordinate of the center of mass of wave function ip. The key observation is that, 
because the additional forces in the accelerated frame of reference depend only on the electron 
coordinates, the additional terms in the transformed action can be expressed explicitely in 
terms of the density. Clearly, these additional terms would have had exactly the same form 
if we had considered the transformation of the non-interacting action functional 

Mp] = l^Mm^^ - f\Mt))dt (12) 

where ipo is the wave function corresponding to the density p(r, t) in a non-interacting 
system, and T is the kinetic energy operator. 

Now, according to Runge and Gross the interacting and non-interacting action function- 
als are related as 

A[p\ = Ao[p\ + Ah[p\ + AM (13) 

which constitutes a definition of the "exchange-correlation" part of the action functional. 

Since the Hartree part of the functional 

is manifestly invariant under the transformation p — > p' we conclude, comparing the trans- 
formation equations for A and Aq that A^c is invariant under this transformation: 
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A.c[p'] = Axc[p] (15) 

This is the most important result in this paper. The tranformation law for the exchange 
correlation potential is derived in the following way. We consider the difference between the 
xc action calculated for two neighboring densities p and p + bp: 

A.c[p + bp] - A,e[p] = f^^ dt' J df'V^.lp- r', i!\bp{^, i!) (16) 

which follows from the definition of V^c as a first functional derivative of A^c with respect 
to the density. We can write the same relation for the transformed densities p' and p' + bp': 

A.c[p' + V] - ^xc[p'] = rft' / f^r-^K,[p'; r-^, ^\bp\^ , t') (17) 

Now, according to eq. (|^) the left hand sides of these two equations must be equal. It 
follows then that the right hand sides are also equal. By doing a change of variable in 
one of the two integrals and using the fact that bp is arbitrary, we easily establish the 
transformation property 

Kc[p'; r, t] = V,,[p- r + x{t),t] (18) 

We now show that an exchange-correlation potential that correctly transforms according 
to equation (|18D automatically satisfies the "harmonic potential theorem" . We suppose that 
for t < to the density of the interacting electron system is described by the usual Kohn-Sham 
equation with a time-independent external harmonic potential Vo(r*), and a static Hartree 
potential VH{po]r) and exchange correlation potential Vxc[po',r\. po{r) is the unperturbed 
ground-state density, given by the sum of the squares of the N lowest lying eigenfunctions 
of the Kohn-Sham equation 

[T + Vo{f^ + VhW, r1 + K.c[po; ^]0.(r) = e,</.,(r). (19) 

The time dependence of these orbitals is given by 0i(r, t) = (j)i{f)exp[—ieit], and they satisfy 
the time-dependent equation 
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[z- - f - Vo{t) - VhW, r1 - Kefpo; t) = 0. (20) 

The system is now perturbed by a uniform time-dependent electric field E{t)\ the center of 
mass moves according to eq. (|l]) and the harmonic potential theorem says that the density 
must evolve according to eq. @. Our task is to show that in fact p{f,t) = poif— Rchiit)) 
is a self-consistent solution of the time-dependent Kohn-Sham equation in the presence of 
the driving electric field: 

[^^ - f - Vo(r) - Vh[p; r, t] - V^,[p- f, t] - e^(t) ■ rl0f (f, t) = (21) 

To prove this fact, we observe that the Kohn-Sham equation, in the presence of the driving 
field, eq. (^), can be generated starting from the Kohn-Sham equation without the driving 
field, eq. (^), simply by applying to the latter a transformation to an accelerated frame 
of reference with x{t) = —RcM{t) where RcAiit) is the solution of the equation of motion 
(H). In fact, using eq. (|35D of the Appendix, it is easy to prove that the transformed wave 
functions (j)[{f^t) = U(j)i{r,t) satisfy the equation 

[^^ - T - Voir- RcM{t)) - Vh[po; r - RcMit)] - Kc[po; r - RcMit)] 

+mdcM{t) ■ f+ -^{vcMit) X B) ■ r + H^i^M^ + ^(5 x RcM{t)) ■ vcM{t)]mr,t) = (22) 

Now, we substitute the harmonic potential form for 

Voir- Rcm) = K)(r) -RcM-K-r + ^Rcm ■ K ■ Rcm, (23) 

and we use the transformation of the Hartree and exchange- correlation potentials Vxdpo] r — 
Rcm) = K.'c[p;^- Using the equation of motion essentially equivalent to the Kohn-Sham 
equation (^) in the presence of the driving field. The only difference is a time depen- 
dent additive term in the potential which can be eliminated by further multiplying each 
wavef unction by the phase factor exp[iS{t)] where 

S{t) = / lirvlMit) - i^RcM{t) ■ K ■ RcM{t) + -A{RcMit)) ■ vcM{t)]dt (24) 
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is the classical action for the motion of the center of mass. Since the new orbitals (with 
the additional phase factor incorporated) yield the density poif — Rcuit)) and satisfy the 
time-dependent Kohn-Sham equation, in the presence of the driving field, with Hartree and 
exchange- correlation potentials evaluated at that same density, we conclude that they are 
the self-consistent solution of the driven Kohn Sham equation. Any approximation for V^c 
that satisfies eq. (|l^) will automatically satisfy the HPT. 

Consider the adiabatic local density approximation [^. The action functional in this 
approximation has the form 

A^ccW, t)] = - J'' dt J dfeMr, t)) (25) 

where e^dp) is the exchange-correlation energy density of the uniform electron gas of density 
p. This expression is manifestly invariant under the transformation p ^ p' since the latter 
amounts to a simple change of variable in the space part of the integral. Therefore, eq. 

On the other hand, consider the linear response theory of Gross and Kohn [|2[ for V^c- It 
has the form 

V;,[p; f, t] = V;^,^^[po(r1; r1 + T /.,[po(r); t - t']Sp{r, t')dt' (26) 

J to 

where it is assumed that the density p{f,t) = po{r) + 6p{r,t) deviates only slightly from 
the static equilibrium density. V!j^^'^[po(r); ^ is the usual local density approximation for 
the static density po. Notice that the kernel f^c is a function of a time difference, but only 
one position (i.e. we have locality in space but not in time). In the case of harmonically 
confined electrons we know from eq. (0) that 5p{f,t) = —Rcuit) ■ Vpo{f) if RcNiif) is 
small. On the other hand, from eq. (pTSD, we know that K..c[p;r, t] = Vxcipoir — RcMit)] ~ 
K'c[Po, r\ - RcAiit) ■ VK-c[Po; ^, and V^^po; r\ = V^J^^[po] r\ . Substituing this in eq. (|26|) 
we see that the integral on the right hand side must equal VV!j.^'^'^[po, r] ■ RcMif)- Because 
this must be true for an arbitrary (small) driven motion of RcAiit) we see that the only 
admissible time dependence of fxc{t " t') is proportional to a 5 function of t — t', i.e. the 
only admissible approximation in this class must be local in time as well as in space. Since 
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the Gross-Kohn potential is non-local in time, it violates eq. (|T8| ) and therefore also the 
harmonic potential theorem. 

How can one make sure that approximate forms of the exchange-correlation potential 
satisfy eq. (|1^ )? A general way is to start from an action that depends only on the "relative 
density" defined as follows. For a given density p we construct the position of the center of 
mass, and then we refer p to an accelerated frame in which the center of mass is at rest. 
The so defined density 

Prei{r,t)= p{f+RcM{t),t) (27) 

is what we call the relative density. By construction, it is invariant with respect to transfor- 
mation to an arbitrary accelerated frame of reference. Therefore, if the action functional is 
written as a functional of the relative density Axc[p{t^ t] = A^clPreii"^, t)] it will automatically 
satisfy the symmetry embodied by eq. ([T5|). The exchange-correlation potential will then 
be given by 

VMr,t] = ^-^l^ir-RcM{t),t). (28) 

OPrel 

An example of application of this method is the procedure proposed by Dobson to go 
beyond the adiabatic linearized LDA without violating the HPT. The idea is to divide the 
density variation into two parts: one corresponding to a rigid shift of the center of mass, the 
other including the remaining internal deformations. The Gross-Kohn frequency dependent 
kernel is applied only to the relative density variation, while the center of mass variation is 
still treated in the ordinary ATLDA. 

In the language of this paper, this idea corresponds to the following procedure. Write 
the relative density as prei{r,t) = po{r) + 6prei{r,t), where po{r) coincides with the absolute 
initial density, because the center of mass is initially at the origin of the absolute coordinates. 
Of course, the correction 6prei vanishes in the case of uniformly driven harmonically confined 
system, but it need not vanish in more general cases (for example, in a non-uniformly driven 
harmonic system). Approximate the action functional as 
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A^J^Prei{r,t)] = - f dt f dre^c{po{r)) 

Jto J 

+ 1 dt r dt' f dfU[po{r); t - t']6prei{r, t)5prei{r, t'), (29) 

where fxc[po,t — t'] is the Gross-Kohn exchange-correlation kernel. The exchange-correlation 
potential is then constructed by taking the functional derivative of d^c with respec to 
Prei{r,t), and translating the position vector from r to r — Rcuif)^ as indicated by eq. 
(p8|). RcMit) itself must be determined self- consistently from the knowledge of the full 
density, according to eq. 

We note, in closing, that all our result could be easily generalized to the case of a uniform 
time-dependent magnetic field. The HPT generalizes to this case. The only difference is the 
appearance of the additional electric field {e/2c)dB(t)/dt x x(t) in eq. (H). 

Acknowledgements I acknowledge support from NSF Grant No. DMR 9100988 and 
hospitality at the ITP where part of this work was done, under NSF Grant No. PHY89- 
04035. 

Appendix - Proof of eqs. (6) and (7) 

We start from the Schroedinger equation eq. (^) for ip and apply to it the unitary trans- 
formation U{t) defined in eq. in the following manner 

f/(t)(z^ -Ho- Vit))U-\t)Uitmt)) = (30) 

We then observe that 

mz-U-\t) = z- - [md{t) + ^^{v{t) X B)] -Y^f, 

^ e - 

+ Y.P^■ + Nmv'it) - -(vit) X B) ■ xit), (31) 
i=i 

and 

U{t)f^-\t)=f, + x{t) (32) 



U{t)piU-\t) =pi + mv{t) - ^B X x{t) (33) 
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The density associated with the wavefunction ip' = U{t)ilj clearly is p'. Furthermore 



U{t)[Ho + V]U-\t) = Ho + Y. yi^i + m,t) + + ■ ^(t) + 7^ ■ i^it) X B) 

(34) 



i=i 2 i=i 1=1 



Therefore, the transformed Schroedinger equation has the form 

[^^ - ^0 - E ^(^* + ^) - "^«(^) ■ E 
<^'' 1=1 i=i 

N 2 
P _ — » Tflll P — » 

-- Er1 • (^T(t) X 5) + + -(5 X x{t)) ■ t7(t)] = (35) 

This can be seen as the time evolution of the wavefunction in the presence of a trans- 
formed time-dependent potential given by eq. (^ ) apart from an additive time dependent 
constant. This completes the proof of eqs. (6) and (7). 
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